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Abstract 

We study the critical behavior at nonzero temperature phase transitions of an effective Hamil- 
tonian derived from lattice QCD in the strong-coupling expansion. Following studies of related 
quantum spin systems that have a similar Hamiltonian, we show that for large N^. and fixed g^N^, 
mean field scaling is not expected, and that the critical region has a finite width at = oo. 
A different behavior rises for A'^^ oo and fixed Nc and g'^/Nf, which we study in two spatial 
dimensions and for Nc = 1- We find that the width of the critical region is suppressed by 
with p = 1/2, and argue that a generalization to A^c > 1 and to three dimensions will change this 
only in detail (e.g. the value of p > 0), but not in principle. We conclude by stating under what 
conditions this suppression is expected, and remark on possible realizations of this phenomenon in 
lattice gauge theories in the continuum. 

PACS numbers: 11.15.Me,11.15.Pg,12.38.Gc,25.75.Nq 
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I. INTRODUCTION AND SUMMARY OF RESULTS 



Large- expansions are useful to approach nonperturbative dynamics of field theories and 
are often used to study their phase transitions at nonzero temperature T. Here we focus 
on the dependence of the critical region of strongly coupled lattice QCD on the number of 
colors, Nf.., and flavors, Nf. This region of temperatures is where mean field (MF) theory 
fails, and fluctuations, that cannot be ignored, drive the system to a nontrivial fixed point. 

What motivates us are the studies of the suppression of this region in the large- Gross- 
Neveu and Yukawa models [1-3]. In particular, the authors in [3], suggest that the same 
phenomenon occurs in QCD at large- A^c, and discuss its physical implications. In apparent 
contrast to that, in the letter [4], which presents numerical simulations of strongly-coupled 
lattice QCD with staggered fermions and large Nc < 48, the authors claim to find that MF 
fails at large-A^c) and as evidence show that the critical region does not shrink with in 
their simulations. 

To understand this apparent contradiction and when to expect a suppression of the critical 
region in general, we study an effective Hamiltonian which is derived in second order of the 
strong-coupling expansion from the lattice QCD Hamiltonian, and that describes the low 
energy effective excitations of mesons. We work in d spatial dimensions and with Nc colors 
and Nf flavors of naive fermions. This gives us flexibility to explore a large parameter space, 
and enables us to show that which large- makes a system trivial depends on which A^ is 
large. Following studies of similar Hamiltonians in the condensed matter literature [5-9], we 
explore the critical behavior of this effective Hamiltonian in three different large- A^ limits: 
(I) the large- A^c limit, where the 't Hooft coupling g'^N^. is kept fixed. (II) the combined limit 
of large- A'c and large- A^j, where a la Veneziano, we fix both N^/Nf and g'^Nc, and (III) the 
limit of large Nf but with fixed Nc and g"^ /Nf. We summarize the results of these studies 
in the next paragraphs. 

In the large- A^c limit we show that the largeness of Nc serves to suppress quantum fluctu- 
ations, and leads to a classical Hamiltonian. This Hamiltonian is a generalized Heisenberg 
antiferromagnet, whose ordered ground state corresponds to the chiral broken phase. A MF 
ansatz in terms of the 'spins' is exact only at T = 0, and is not adequate to study the 
critical behavior near Tc- In particular, the critical exponents are determined by Nf and 
d, and are not expected to belong to the Gaussian fixed point. One can, as a first step, 
analyze the transition in MF, but by definition, this analysis is bound to fail within the 
critical region, which we show to have a nonzero width at A'c = oo. Unfortunately, since at 
Nc — oo the temperature Tc turns out to be infinite, then it is difficult to learn about the 
behavior of transitions in planar QCD from the 't Hooft limit of our effective Hamiltonian, 
and we therefore use these results only to understand [4], but discuss how they should be 
improved. 

In the combined limit of large Nc and Nf, the situation is different, and one can solve 
for a MF ansatz that is exact at Nf, Nc ^ oo for all T [Tc is finite is this hmit). Reviewing 
known results, we show that this solution yields a correlation length that diverges with a 
MF exponent only for > 4, and otherwise for 2 < c? < 4. Extending this, we also show 
that the width of the critical region depends on Nf, and Nc only through their ratio Nc/Nf, 
and is nonzero. 

In the limit of large- A^/ and fixed N^ a MF ground state is again exact for all T, but 
leaves chiral symmetry intact even at T = 0. In d = 2, a global minimum of the system is the 
"spin-Peierls" state, which breaks lattice translations and rotations [7, 8]. A corresponding 
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analysis ind = 3 includes an extensive search in the space of all possible ansatzc for the global 
ground state, and is out of the scope of this work. Since this limit has a restricted relevance 
for QCD (for example, asymptotic freedom is lost here) but is crucial for our purpose (see 
below), we proceed to study the restoration of these lattice symmetries at finite T for the 
d = 2 system, and for simplicity also fix iVc = 1. We believe that a generalization to d ~ 3, 
and A^^c > 1 will change the results only in detail but not in principle, and find no point to 
do so. (For example, see the Nc> 1 generalization at T = 0, and d = 2 in [7]) 

The lattice symmetries are restored in a second order transition at the finite Tc, whose 
critical behavior is closely related to the behavior seen in [1, 3, 10, 11] for the systems studied 
there. We show that taking Nf — > cx) before t = \T — Tc\/Tc 0, makes MF exponents 
exact. Switching the order of the limits, we find that the Landau-Ginzburg- Wilson (LGW) 
action for the transition describes scalar fields coupled through 0{1/Nf) interactions. As 
emphasized in [3, 10], this suggests a crossover behavior, which we find to occur when 

t ~ 1///V/. 




Combining the results of [3. 10. 11], and of our study here, we emphasize that the critical 
region in large-A^ phase transitions is suppressed only when the LGW effective action for 
the transition has an overall factor of N"", a > 0. This suppresses thermal fiuctuations in 
the order parameter, and can make MF scaling exact. Our message here is that this does 
not happen in all large- A/" treatments, and in particular does not happen for the chiral phase 
transitions of the effective Hamiltonian that we study here. 

We begin this report in Section II where we introduce the effective Hamiltonian of lattice 
QCD in the strong-coupling limit. We then move to discuss large and fixed Nf in 
Section III, large Nc and large Nf in Section IV, and the large Nf with fixed Nc in Section V. 
We conclude in Section VI, and make proposals for future research in Section VII. 

II. THE STRONG-COUPLING LIMIT: THE EFFECTIVE HAMILTONIAN 

Strong coupling expansions were used in lattice gauge theory since its early days. In 
particular, they were performed for the gauge group SU{Nc) as well as for SU{3). For 
example, in the pure SU{Nc) gauge theory, expansions were made for the string tension 
and the for free energy as well as for the beta function [12]. In the case of QCD with Nc 
colors, and various types of fermions, calculations were made for the effective low energy 
action of hadrons and for mesons masses [13, 14]. The result of these treatments (which is 
most relevant to the discussions is Sections III-IV) is that for large values of N^ the natural 
expansion parameter in the strong coupling series is the inverse 't Hooft coupling l/[g^Nc) 
rather than (see also [15, 16]). More precisely, in these expansions one fixes A^^^ and take 
g'^Nc 3> 1 to find that quantities which are 0{N^) at large- iVc, such as mmeson, ''^baryon/A'c, 
fl/Nci etc, can be written as a power series in {g'^Nc)~^ with coefficients that depend on Nc 
and that become finite at A^c = oo- 

The result of this procedure may be different from taking the opposite order of the limits 
(first A'c — > cxo and only then g'^Nc — >■ oo). As discussed in [16, 17], one may find that these 
limits commute only as long as {g'^N^Q < g'^Nc < oo. Indeed for the two-dimensional pure 
gauge theory with the Wilson action one finds {g'^Nc)o = 2 [17]. For higher dimensions the 
situation is more involved but there still is a finite range of g'^Nc in which the theory is in a 
strong coupling phase (For example see early analytical results in [18] and recent numerical 
results in [19]). Since fermions are quenched in the 't Hooft limit, we take these results to 
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apply for QCD (with a fixed nonzero mass) as well, but we are not aware of any analogous 
results in the Veneziano limit of fixed Nf/Nc- 

In this paper we follow the former approach (which is also the approach of Euclidean 
treatments such as [4]): we fix Nc and use an expansion in [g'^N^'^ <^ 1. The result of this 
expansion in the Hamiltonian formalism is that at second order one obtains the following 
effective Hamiltonian, 

J 



n+A' 



(2.1) 



which is derived for naive fermions in [13, 20]. Here n denotes lattice sites, fx = Xi, . . . 
denotes the lattice directions, and the couphng J scales with the 't Hooft couphng like 
(g^Nc)'^ [13, 21]. The operators generate the U{ANf) algebra 



(2.2) 



and are defined as 



(2.3) 



af3,a 



Here the fermion indices a and /? are Dirac-flavor indices that range from 1 to 4iVj, while a 
is a color index that ranges from 1 to Nc- The matrices are the generators of U{ANf) 
in the fundamental representation, and obey 



trM''M'' = ^V' 



'Pi- 



(2.4) 



The Hilbert space on site n is an irreducible representation of U{ANf) that corresponds to 
a rectangular Young tableau with N,, columns, and m rows, see Fig. 1. The corresponding 



m 



FIG. 1: The representation of U{4:Nf) carried by Q^. m is related to the baryon number at the 
site according to m = B + 2Nf, with \B\ < 2Nf. 



baryon number on such a site is -B = m — 2Nf. In this work we restrict to configurations with 
zero average baryon number, and in particular examine configurations with baryon number 
B on the even sites and —B on the odd sites. This means that we put conjugate U{ANf) 
representations on adjacent sites, represented by Young tableaux with m rows on the even 
sites, and with (A^ — m) on the odd sites. Following [5-7] we take m = 2Nf in Sections III, 
and V, and m = 1 in Section IV. The reason is that these cases are the simplest to analyze. 

The Hamiltonian Eq. (2.1) is a generalized Heisenberg quantum antiferromagnet, and 
its Neel ordered state spontaneously breaks the global U{ANf) symmetry. Adding explicit 
symmetry breaking terms can reduce U{ANf) to its subgroup SU{Nf)L x SU{Nf)R xU{1)b 
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(e.g. see [20, 22]), and as a result, the Neel phase corresponds to the spontaneous breakdown 
of chiral symmetry.^ 

Note that Eq. (2.1) is a result of Rayleigh-Schro dinger degenerate perturbation theory 
applied to split the degeneracy of the zero gauge-flux sector of the strong-coupling lattice 
QCD spectrum. Excitations of strings with a nonzero gauge flux cost an energy that scales 
like ~ g'^Nc ^ 1 and are neglected here. As a result the chiral phase transition of Eq. (2.1) 
will not involve a condensation of Polyakov lines of the type discussed in [23]. Adding the 
effects of these gauge strings will make the description of the transition more realistic, and 
in particular will relate chiral symmetry restoration to deconfinement. This restriction to 
the zero flux sector means that the Hamiltonian Eq. (2.1) describes chiral dynamics of QCD 
only at temperatures T < g'^Nc- As a result one should be careful when 'extrapolating' 
the behavior of Eq. (2.1) (or its Euclidean counterpart) to QCD at temperatures close to 
the transition. This is a significant difficulty that exists in the strong coupling approach to 
finite temperature lattice QCD, which we cannot resolve here, and that the reader should 
be aware of. 

Nonetheless, since the model in Eq. (2.1) shows different critical behaviors in different 
large- limits, it is a sufficient and useful choice for our study, and can tell what conditions 
are needed for a large- A^ transition to have MF critical exponents. Also, the Hamiltonian 
Eq. (2.1) is sufficient to understand the work of [4], which also neglects excitations of the 
gauge-flux. 



III. THE LARGE-AT^ LIMIT : A SIGMA MODEL 

As mentioned in the previous section, we now fix the 't Hooft coupling g^Nc (and hence 
J) and take the large-A'c limit of Eq. (2.1). Writing the partition function of Eq. (2.1) 
using generalized spin coherent states gives a non-linear sigma model (NLSM) that was first 
derived in [7]. The a field at site n is a 4Nf x 4Nf hermitian, unitary matrix given by the 
U{4:Nf ) rotation, cr„ = UnAU^, of the reference matrix A. In this section we work with 
m — 2Nf that fixes 

(3.1) 



A = 



^2Nf 



-l2N, 



and makes (7„ an element of U{4:Nf)/[U{2Nf) x U{2Nf)]. The action of the NLSM is 



2 Jo 



dt 



Y;^TrAUidtUr, + 2x^Tr ( 



(3.2) 



where x = J /AT gives the coupling in units of the temperature T. Finally note that the 
global U{4:Nf) transformations are realized here as Un — > VUn, or (T„ — > V^cr^V^'^. 

For large values of Ac, and x ~ 0(1), the overall factor of Nc/2 suppresses fiuctuations, 
and the MF anstaz 

/ \ _ r +^ ^ ^ even, , . 

\'^r.) - [ -A ne odd, ^"^-"^^ 

which breaks U{4:Nf) to U{2Nf) x U{2Nf), is exact. An expansion around the MF ansatz, 
results in a chiral theory of mesons with 0(1/ N^) interactions [13]. 



^ See, however, the remark at the top of Section IV. 
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In contrast, if x is of 0{1/Nc) then fluctuations in the second, 'interaction' term, are not 
suppressed, and many sigma field configurations contribute to the path integral. Nonetheless, 
the largeness of Nc suppresses quantum fluctuations for any x, by making dfU ~ in the 
first, 'kinetic' term. This can also be seen by rescahng the operators — > Nf. Q'^ in Eq. (2.2), 
to find that the commutation relations vanish at N^. — oo. In fact, this is in analogy with 
the large-^ limit that makes quantum antiferromagnets into classical systems. In our case 
we get the classical Hamiltonian 

"^classical .rp ^ {(^n^^n+fl 

whose critical temperature Tc ~ 0{NcJ). This means that x scales hke ]^ x and therefore 
that at T ~ 0{T~), x is of 0{l/N^. As a result the MF ansatz Eq. (3.3) is not predictive 
at these temperatures, and treatments such as [27], are expected to fails there. This can be 
also seen in the following way. The partition function of this classical system looks like 



Z = Da exp 



Tc 



(3.5) 



where a is a pure number of 0(1) given by NcJ/{4:Tc).^ In terms of T/T^, Eq. (3.5) has no 
Nc dependence at all, which means that the critical region has a finite width at Nc — oo. 
Only outside this region will MF be a good description, like it is for ordinary spin systems. 
MF theory does not fail here since it is not supposed to work at all, and we believe that this 
is what [4] saw numerically. One might think that MF scaling fails because the ansatz used 
above is temperature independent, and as emphasized above, is strictly exact only for low 
T. In the next section we show that this expectation is too naive. 

Finally we note that the estimate of Tc ~ Nc has been obtained in other strong-coupling 
expansions as well (for example see [4, 24]). At first sight this is puzzling since there 
should be only one scale in this theory which is IVqcd- This situation is reminiscent of the 
one discussed in [25], where the author shows that the thermal corrections to {^pip) in the 
framework of continuum chiral perturbation theory are of 0{T'^/Nc). Naively, this leads to 
expect Tc ~ ^/Nc- This 'paradox' is resolved in [25] by noting that the expansion in T must 
have a finite radius of convergence of the order of the Hagedorn temperature. The latter is 
an external notion to the chiral lagrangian, and can be thought of as proliferation of string 
of nonzero gauge-flux. 

Indeed, the action of our nonlinear sigma model Eq. (3.2) is the chiral lagrangian in our 
context, and is reached only when one neglects fiux excitations (see discussion in Section II). 
As a consequence, it seems that the reason that in strong-coupling Tc ~ Nc, and not a 
function of just g'^Nc, is the discard of flux condensation. In fact the authors in [26] study 
the interplay of deconfinement and chiral symmetry restoration in strong coupling, and find 
that an assuming fiux condensation gives a chiral symmetry restoration temperature which 
is 0{N^) at large- A^^c- We thus expect that the inclusion fiux condensation will lead to a 
Tc ~ 0{N^) in the Hamiltonian approach as well, and that it will be the critical behavior of 
the resulting system which may be indicative to what happens in the 't Hooft limit in the 
continuum. 



^ For example, in [4] one finds that the temperature is given by T^/Nc = 1.5525(3) -|- 0{1/Nc). 
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IV. THE COMBINED LARGE-AT^ AND LARGE-A^/ LIMIT : SCHWINGER 
BOSONS 



The aim of this section is to investigate a case where the MF ansatz is T-dependent 
and exact at A^" = oo. This can be achieved by taking both Nf, and Nc to be large, while 
keeping their ratio fixed {g'^Nf. is still kept fixed here). This limit was first studied by [6], 
for m = 1 with the Schwingcr bosons method, and leads to a spontaneous breakdown of 
U {4:Nf ) U{1) xU {4:Nf — 1). This tells us that chiral symmetry may be realized differently 
than in QCD, and that to break U{ANf) in the same pattern as in Section III, one needs 
to generalize the procedures in [6] to m = 2A^/. This was done in [9], which found that the 
simplest ansatz for general m gives the same MF equations as for m = 1. This encourages 
us to choose m = 1 as well, and to postpone a discussion of a more sophisticated ansatz to 
possible future research. 

We now move to review the Schwinger bosons method, and present some of its results. 
For a more detailed discussion, we refer to [6, 28, 29] and, in the context of QCD, to [21]. 
The first step is to write 

or? _ / ^« ■ ■ n e even, , . 

^""^ \ hl-{-M^)* -h^ neodd, ^^■^> 

where the operators h^n are bosonic fields that live on the lattice sites, and have only Dirac- 
flavor indices. This is an acceptable representation of Q'^, since ban obey ban, b^rn ~ 
Sap(^nm, ^ud therefore Eq. (2.2) is respected. To set the representation of the operators in 
Eq. (4.1) one puts Nc bosons on each site 

ANf 

Y: bijar. = Nc. (4.2) 

This constraint, together with the fact that the bosonic single-site wave function is symmetric 
in Dirac-fiavor indices, sets the U{4Nf ) representation of all lattice sites to be the tableau 
of Fig. 1 with m — 1. Using Eq. (4.1), one can now represent the partition function, Z, with 
bosonic coherent states, and obtain the path integral 

Z = J' DbDb* ex.p{-A), (4.3) 



rl/T 



Yl Kn^rban + KnKn+fibpnb)3n+ii 



(4.4) 



with J = JNf / {2Nc). Since we keep Nc/Nf and J fixed, then J is fixed as well, and we 
drop the bar henceforth. 

The prime in Eq. (4.3) means that the path integral is constrained to obey Eq. (4.2). 
Adding a Lagrange multiplier, A^, to keep this constraint, and a Hubbard- Stratonovich 
(HS) link field, Qrip,, to decouple the quartic interaction, one finds 



/ 



dr 



iV f 



(4.5) 
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An integration over the bosons, gives the effective action for Q and A, 



+ trlogG-^L (4.6) 



where k = Nc/ (4iV/), and where G is the propagator of a single boson, such that log det G ~ 
0(1). For large values of Nf, but k, ~ 0(1), and T/J ~ 0(1), fluctuations are suppressed 
around any stable MF ansatz. The action of the ansatz Qnnir) — Q,\n — i\ is 



IMF 



ANfNs 
T 



4J V 27 Ns(5 



(4.7) 



where j3 = 1/T. Also, Ng is the number of sites, k takes values in the Brillouin zone 
of the even sublattice. Defining 7 = ^J2i=iCos{ki/2), one finds that the function ujk — 



— AcPQ'^^l becomes a;| ~ + c^fc^ at low k, with the mass gap A = ^/A^ — Ad?Q'^, 
and with = dQ"^. A minimization of Eq. (4.7) with respect to A and Q, gives the MF 
equations 



1 2 ^ A / . . 1^ 



dJ fe ojk 

with nsiuj) = {e^"^ — 1)~^. These equations are exact at large- A", and their solution gives the 
large- A^ phase diagram. This is different from Section III where the mean-field ansatz is exact 
only at low T. In fact, we believe that the Schwinger bosons approach is an appropriate 
candidate to be the 'finite T mean field theory' mentioned in [4]. Another candidate is 
presented in Section V. 

The point at which the inverse correlation length A vanishes, signals the condensation of 
the bosons, which breaks U{4Nf) to f/(l) x U{4Nf — l). Investigating Eqs. (4.8)-(4.9) shows 
that this happens for T = at k > Kc,^ and that the symmetry is restored at T = Tc{k), 
which is finite [6, 29]. 

To extract the critical exponents, one expand Eqs. (4.8)-(4.9) around T^. We illustrate 
this in Appendix A, where we redo the calculation of [29] to show that MF scahng holds 
only for d > 4, and that in terms of the reduced temperature t = \T — Tc\/Tc, one finds that 
A ~ t^/^. For 2 < d < 4, infrared (IR) fiuctuations change this behavior to t ~ — A^logA 
for d = 4, and to A ~ t^/i^-"^) if 2 < d < 4, which coincides with the behavior of the CP/v 
model at large- A^ [30]. Extending these results, we also show that the critical region for 
2 < c? < 4 is nonzero and depends on Nf and N^. only through k. 

We conclude this section by noting that the overall factor of 4A^/ in Eq. (4.6) does not 
contradict the fact that IR modes are not suppressed. The simple reason is that the fields Q 
and A are U (4A/) singlets, and are not the order parameters of the transition. This means 
that Eq. (4.6) is not the LGW action and that fluctuations in the order parameters need not 
be suppressed. Large-A^/ suppresses fluctuations in Q and A, which means that the results 
one obtained here are exact ai Nf — 00. Taking these fluctuations into account leads to 
corrections of 0(1/ A^/) to the critical exponents, and to Tc- 



3 Kc is ~ 0.19 for d = 2 [6], and ~ 0.0778 for d = 3 [28]. 
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V. THE LARGE-A^/ LIMIT: FERMIONS 



Here we take Nc = 
to write Eq. (2.1) as 



1 and work in two spatial dimensions. We begin by using Eq. (2.4) 



a/3 



J 



/ na 



(5.1) 



where we define J = JNf/2 ~ Nf/g^. Since we take large- A^/ with fixed g"^ /Nf then J is 
fixed as well, and for brevity, we drop the bar henceforth. To set the representation of the 
operators in Eq. (2.3) we follow [6, 9] and put 2A^/ fermions on each site 



2Nf. 



(5.2) 



This also means that the second term in Eq. (5.1) is a constant which we drop. Since the 
wave function of the fermions is an^i-symmetric in the Dirac-fiavor indices, then Eq. (5.2) 
puts it in the representation given in Fig. 1 with m — 2Nf. The partition function of 
Eq. (5.1) is [7] 

/ttT r _ 

^DXJ DQDQ*DipDiPe-^, (5.3) 

K n n,iJ. "-^ n,n n ) 

(5.4) 

As in the previous section, A„ is a constraint field that keeps Eq. (5.2) on each site, and 



Qnfj.{T) is a HS field that decouples the four- Fermi interaction of Eq. (2.1). Apart from 
the global U{4Nf) symmetry, Eq. (5.4) is also invariant under the U{1) gauge transfor- 
mation ^ e^^"(^Vn,^r. ^ i'r.e-'^-^^\ Qr^^. ^ e'^^^^^^'^^+f^^^y^ Qr^f,, \r^ ^ A„ + ^^A^. 

Demanding that the fields Q, will be single- valued, and that A„ e [— ttT, ttT] remains 
time-independent, restricts A„(r) to be time-independent as well. 
To proceed, one integrates over the fermions and obtains the action 



Aff = 4Ar 




77 E \Qnfii.\^ + 9 E -^^ 



+ tr logL'-^(g,A) 



(5.5) 



where D{Q, A) is the Dirac operator of a single fermion, such that tr logD ~ 0(1), and can 
be read from Eq. (5.4). For T/Jr^ 0(1), the overall factor of 4A^/ in Eq. (5.5) suppresses 
fiuctuations in Q and A, and one can exactly solve for the ground state which was shown 
to be the spin-Peierls state [5, 7, 8]. It has A„ = on all sites, and \Qnfj.iT)\ = q on a. 
subset of the lattice links B, but zero otherwise. This ground state is four-fold degenerate, 
and is shown pictorially in Fig. 2. It is clear that this state breaks lattice translations and 
rotations, which are then restored at the critical temperature Tc. To discuss this transition, 
one writes the MF effective action for the spin-Peierls state 



^MF — 7?, 



11 1 

—q^ X - X 2 + —lognF{q)nF{-q) 



(5.6) 



9 



A B 



: [ [ ] ] 
■ [ ] ] ] 

; [ ] ] ] 












: [ ] ] ] 
: [ ] ] ] 







C D 



FIG. 2: The four degenerate spin-Peierls ground states, which break translation and rotation 
symmetries. The low energy effective action we compute begins from state "A". 



where /3 = 1/T, and npie) = (e^^ + 1) ^. A minimization of Eq. (5.6) with respect to q 
yields the MF equation 

g/2J = tanh(g/2r), (5.7) 

that describes a stable ordered phase with q > below Tc = J, and a MF scaling q ~ 
(Tc — T)^/^ in its vicinity. Here we reach the same 'puzzle' as did [1, 3, 10, 11] for the Gross- 
Neveu, Yukawa and classical 2D Heisenberg models: the result above is exact at Nf = oo, 
but for any finite Nf it is wrong. There, the scaling of q is dictated by the symmetry 
breakdown, and by d. The resolution of this puzzle is that the critical region shrinks with 
Nf [3, 10, 11]. We now turn to show how this happens for our model. 

Wc begin by calculating the effective action for fluctuations by replacing \n A„, 

and Qnjj, q+ Qnfi in ^eff (Here q is the solution of the MF equation, and is nonzero 

only on B). We expand Eq. (5.5) to 0{1/Nf) in Appendix B, where we work in Matzubara 
space and show that the masses m of the zero Matzubara flelds (f) = Q{uj = 0) vanish at Tc 
like ~ t = |Tc — T\/T. All other fields are massive, and we proceed by integrating them 
out in Appendix C. The effective action we find is 



*/ {n},fiu 



J {n},ij.iy 



(5.8) 



with m^,f ~ 0(t), V^^^ ~ 0(t^/^), and V^'^^ ~ 0(1) (detailed expressions for these are given 
in Appendix C.) 

As emphasized by [3, 10, 11] it is now clear where does the puzzle come from. Taking 
A^/ — > oo before i — > gives a Gaussian model. Switching the order of the hmits, we get a 
weakly coupled (f)^ theory, and the four-fold degeneracy of the spin-Pcicrls state leads one to 
expect that the universality class of the transition is of a Z4 model in c? = 2 (see for example 



10 



[31] and its references). This points to a crossover behavior, where the susceptibihty x 
diverges hke 

x = tN';, (5.9) 

and where f{x) is a scahng function that determines the critical behavior. 

A calculation of f{x) is simpler in the high temperature phase, where q = 0. There, the 
effective action has only diagonal and degenerate mass terms, and no cubic interactions. In 
momentum space the action reads 

1 1 

^eff = E l^^^M + E E [^uw({fc})</'fel^</'te,■</•fe3^■</'fe4^ + h.c] , 

ki ^ felfe2 »3 

h^h^ kl 

(5.10) 

where the momentum k belongs to the first Brillouin zone of the even sublattice, and the 
index i — 1,2,3,4 denotes the four links outgoing from each even site, see Fig. 3. In 

• X • X • 

2 

X • 3 X 1 • X 

2 4 2 

• 3 X 1 • 3 X 1 • 

4 2 4 
X • 3 X 1 • X 

4 

• X • X • 

FIG. 3: The lattice of even sites, denoted by X, and odd sites, denoted by •. Each even site has 
four link fields denoted by (pi, with i = 1, 2, 3, 4. 

Appendix C we show that rri^ ^ {T - Tc)/{4TTc), and give the form of F ~ 0(1). 

We calculate the susceptibility Xij{k,k') — {(f)i{k)(l)j{k')) to leading order in 1/A^/ in 

Appendix D, and show that it coincides with Eq. (5.9) if a; = t \f^\[^fi ^'^^ f{^) = 1 — 1/a;^- 
This tells us that only if a; = oo then one obtains the MF scaling of x ~ t~^. As x decreases 
from infinity, this behavior changes, and will eventually be determined by the nearest fixed 
point of the (/f)^ theory of Eq. (5.10).^ This crossover behavior will occur when x ~ 0(1) or 

To make a connection with the results of [3], we recall that a 0^ theory in d dimensions 

2 

and a 0{1/N) coupling enters its critical region when its mass squared is < (l/N)^. 
Our result fits here if we put d = 0, and gives the largest critical region possible. The 
reason is that here the bare propagator of |0| is momentum independent, and is simply 
given by (|m^)^^. This degeneracy will be removed by higher orders in 1/iV/, and can be 
avoided by analyzing other ground states, such as the flux phase [8]. This also means that 
in a generalization to other values of Nc and d, it is reasonable to expect that the width of 
the critical region will be given by i ~ l/N^, with p > 1/2, or by a possible logarithmic 
modification of that. 



Or to a possible modification of that by topological effects, induced by the compactness of the gauge field. 
These have been shown to be important in the T = phase transition between the spin-Peierls and Neel 
states [31]. 
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In contrast to the combined limit of large- iVc and large- iV/, the critical region is suppressed 
here because the order parameter of the spin-Peierls state is the Q field, whose action has 
an overall factor of ^Nj. The latter suppresses the interaction terms in the LGW action 
given in Eq. (5.10), and can make MF scahng exact. 

VI. SUMMARY 

We work with an effective Hamiltonian derived from the lattice QCD Hamiltonian in the 
strong coupling expansion at second order. The large- A^^c limit of Eq. (2.1) is the analog 
of the large-^* limit of quantum antiferromagnets. For the latter, quantum fluctuations are 
suppressed, and one is left with a classical magnet, which is by no means a simple object. In 
particular, its critical behavior cannot be approached with a MF spin ansatz, because the 
latter is a good approximation only at sufficiently low temperatures. The spin S determines 
the energy scale of the system, and for our case of Eq. (2.1), this leads to Tc ~ Nc-, and to 
the fact that in terms of T/Tc, the critical region of our model is finite Nc — oo. We 
suggest that this is what [4] saw numerically in the action formahsm, and that MF does not 
fail there, since it was not supposed to work in the first place. 

The combined limit of large Nc and\aigc-N f, but fixed ratio N,./Nf. is more convenient to 
study the transition. There, the critical temperature is finite, and the MF ansatz, which is 
exact ai Nf — Nc — oo, is temperature dependent, and describes the critical region exactly. 
For c? > 4, it gives Gaussian scaling, but for d < 4, IR modes change this and one can 
calculate the exponents exactly. Also, the width of the critical region depends on Nf and 
Nc only through Nc/Nf, and is nonzero. 

In the large-Nf limit, we have studied the case oi d = 2, and Nc = I. Here the MF ansatz 
is also exact for all T, and breaks lattice symmetries. This transition is characterized by MF 
critical exponents if Nf is sent to infinity before sending t = \T — Td/Tc to zero. In contrast. 



one can take A^^ oo, and t with fixed x ~ tJNf. This leads to a crossover from a 



Gaussian fixed point (at x = oo) to a nontrivial fixed point of a weakly coupled scalar field 
theory. A schematic view of a possible renormalization group flow for this case is given in 



To conclude, it is clear that the critical region in large- A?^ phase transitions is suppressed 
when the Landau-Ginzburg- Wilson action has an overall factor of A^" with a > 0.^ Taking 
the A^ — > oo can then remove IR fluctuations in the order parameter, and makes Gaussian 
scahng exact. We stress that this does not happen for all large-A?" theories, and for our 
model Hamiltonian of Eq. (2.1), it happens only in the limit of large- A^/ and fixed Nc- 

VII. FUTURE PROSPECTS 

A generalization of the large- A^j discussion to = 3, and A'c > 1 is straight forward, 
but relatively costly. For example, one has to find the true vacuum out of all possible 
ansatze (although an effective action for fluctuation around a metastable vacuum can also 
be calculated) . We beheve that the results of such a generahzation will differ from the case 



^ a = 1 in conventional cases where N is related to a global symmetry, while for deconfinement in a pure 
SU {N) gauge theory one expects a = 2. 




Fig. 4. 



12 




FIG. 4: A cartoon of a possible renormalization group (RG) flow in the {t,l/Nf) space {t is the 
reduced temperature — •) RG drives the system away from the Gaussian fixed point denoted 
by the This happens along the axis towards a possible nontrivial fixed point denoted by 

the 'X' (which for simplicity we placed at Nf = 1). The critical region is to the left of the dashed 
line, 1/Nf ~ 

studied in Section V only in detail, and not in principle.^ In view of that, and the limited 
relevance of the large- A^j limit for QCD in our context, we find no point to do so analytically. 
A numerical check of our predictions, and of the possible generalization mentioned above, 
can be made with similar methods to the one of [4]. Ideally this should be complemented 
with an analytical study in the action formalism at large- A^/, in similar lines to Section V. 

As suggested by [25] and [3] the critical region may be suppressed in planar QCD (large- 
Nc with fixed g'^N^ and fixed Nf). While [3] suggests that the outside the critical region, 
one observes a MF behavior, then [25] expects that the chiral condensate will not change as 
a function of T at all until the critical region is reached. It will be interesting to check these 
theoretical expectations, but unfortunately the current study of our model Hamiltonian 
cannot (and is not intended to) accomplish this. The reason is two-fold; First we are 
very far from the continuum limit. Second, the Hamiltonian Eq. (2.1) is not designed to 
explore temperatures which are close to the deconfinement temperature. While the latter 
problem can be in principle tackled by including gauge-flux condensation (see the discussion 
at the end of Section III), it is hard to see how to overcome the former. Approaching 
continuum physics with numerical simulations at large- A^'c and/or large- iVj is very costly. 
While simulating dynamical quarks at large- A^c is unrealistic, then quenching them or using 
dynamical fermions at large- A^/ is expensive. 

This leads us to seek phase transitions in other related models in which the critical region 
may be suppressed as well. An obvious choice, which is numerically 'cheap', is deconfinement 
phase transitions in large-A'c, where the role of the LGW action is played by a Polyakov 
loop action, in the spirit of [32]. The fact that the pure gauge theory has a first order 
deconfinement transition at large- A'c (e.g. see [33]) means that to study such a phenomenon 
in pure gauge, one has to approach a hidden Hagedorn transition, where the mass of the loop 



^ For example, one expects that both x and f{x), depend on t and Nf, in a way that changes with the 
dimension d. 
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vanishes. This was done in [34], but the results are not accurate enough to unambiguously 
determine the critical exponents/ 

Adding bosonic matter to the pure gauge theory may change the deconfinement transition 
into second order. In particular, using bosons in a bi-fundamental representation with Nf. 
even, leaves a Z2 subgroup of the Zj^ center intact. As a result, the deconfinement transition 
(if second order) is in the Ising model universality class, and by adding a fiavor structure to 
the bosons, one can study the chiral symmetry restoration as well [35]. Also, these scalar 
theories are interesting on their own right, as they may serve to check the l&xge- Nc relation 
between supersymmetric and non-supersymmetric gauge theories discussed in [36] . 

APPENDIX A: THE CRITICAL EXPONENTS AND THE WIDTH OF THE 
CRITICAL REGION IN THE COMBINED LARGE AND LARGE iV/ LIMIT 

To extract the critical exponents from Eqs. (4.8)-(4.9), we restrict to the disordered 
phase, where there are no Bose condensates, and where one can replace the momentum 
sums in Eqs. (4.8)-(4.9) with the integrals 

r (dkYj^ coth(?/ VI - 7' + 7^ A^) 
~ J [2^) VI - 7' + 7'A^ ' 

(Al) 

^^^^ ^ J [2^) Vl-7^ + 7^A^ • ^ ^ 

Here we also divided Eq. (4.8) by Q/X, and defined y = /?A/2, and A = 8\/ Jd. Let us now 
concentrate on Eq. (A2), and on values of n where the T = ground state is ordered. For 
each such k, there exists y = ydi^), where Eq. (A2) gives A = 0, signaling the transition. 

Following [29] we denote the integral on the right hand side of Eq. (A2) by I{y, A), and 
write Eq. (A2) as 



7(y,0)-7(ye,0) 




coth(|/ VI — 7^) 

(A3) 



coth(|/ VI - 7^ + 7^A^) 



VI - 7^ + 7^A^ 



For /' = —dl{y, 0)/dy we find 




^ > 0, (A4) 



sinh y^/Y^^J^ 

which means that for y ~ y^, the left hand side scales like {yc — y) > 0. Next, we expand 
the right hand side of Eq. (A3) in A, and note that sufficiently close to Tc the integral is 



The reason is the high calculational cost, and the fact that the field configurations used there are meta- 
stable rather than stable. 
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controlled by IR momenta. For these we write (1 — 7^) = k^/Ad to find that the IR part of 
the integral is 

(4(iA)2 fRfdkV 1 



y Jo \2nJ k^{k^ + AdA^)' 

Here R separates between the IR modes and rest, and we take 1 ^ 2A^(i. As a result 
Eq. (A3) becomes 

d>4, 



iVc - y)/yc = Ad 
with the finite given by 



log^ 



. d-2 



2 < d < 4, 



{Adfr 



A, 



Vc 



1 



'dkYl_ 
271 ) fc4 



R/^ dx 



2%) x^{x^ + Ad) 



ci>4, 
d = 4, 

2 < d < 4, 



We proceed to write Eq. (Al) as 

A(l-A2) = 2« + l-|(^^^ v/l-7' + AVx 



coth(y^l-72 + AV), 

which at T = Tc(«;) gives 

Ac = 2«: + 1 - |(^^ j yr^coth(y, ^W). 



This means that 



with 



A, - A = --A\yJ"' + K) + {y, - y)I", 



I 
I'" 



(dkV 



\2n J sinh^ {VcV^ ~ 7^) ' 



Together with T = AJd/{16y), this gives 



t={T-T,)/T, = BaA^ + A'aX 



R 



d > 4, 



A^log^ d = 4, 
A'^-^ 2<d<4, 



(A6) 



(A7) 



(A8) 

(A9) 
(AlO) 

(All) 

(A12) 



(A13) 
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with 



K = -4,(1 - (PAS) 

We choose to work with vahics of k,, where both A'^, and are positive.^ Also note the 
should include other 0(1) coefficients that come from momenta > E? in the left hand 
side of Eq. (A3). 

For sufficiently small A, Eq. (A13) leads to A ~ t^/'^ for d > 4, while for d = 4 and 

2 < d < A, one finds that —t ~ A^logA, and A ~ t^/i'^-'^)^ respectively. The width of the 
critical region is given by comparing the two terms in the right hand side of Eq. (A13), and 
we find that MP behavior will fail when 



4 



for (i = 4, and when 



l/(4-d) 



for 2 < d < 4. Since B^/A'^ depends only on k we conclude that the critical region of the 
phase transition in the large Nc and large Nf limit does not shrink with Nf or Nc, and is 
nonzero (Here we ignore scenarios where one tunes n to special values where Ba/A'^ ^ 1, 
and concentrate in generic values of k where B^/A'^ 0{1).) 



APPENDIX B: CALCULATION OF A^s IN THE LARGE-ATy LIMIT 



In this section we calculate the large- A^/ effective action for the auxiliary fields Q and 
A around Tc, where the spin-Peierls state dissolves. We begin by replacing Qnp. q + 



-Qnu, and A 



=A„ in Eq. (5.5), and move to Matzubara space, where we denote 



^y^^nfj,: ^^^^ 

the frequencies of the fermions by e, and the external bosonic frequencies by cu. The action 
we get is 



A ^ A 



MF 



(Bl) 



exp(— 5A) = (exp < 



Y^^riYj V'n(e)V'n(e) E ["^"(c + ^)Qniii(^)'il^n+ilie) + h.( 



^4A^ 



^ Numerically, we find that for k ^ 0.21 then is negative for d = 3. This leaves us with the window 
0.0778 < K < 0.21 to explore the symmetry breakdown. We are not aware of any discussions in the 
literature with regards to the cases where A'^ <0. 
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The (, )q in Eq. (B2) means an average with respect to the free fermionic action given by 
Ap = ipDoip, with Dq connecting fermions that reside on edges of a hnk on B, and is given 

by 

= ^ ( ' ] ■ (B3) 

The calculation of 6A is straightforward and the result is given pictorially in Fig. 5, 
where the external legs represent the fields Q and A. We denote the contributions to SA as 




FIG. 5: Summing the bubble diagrams. Internal fermion loops are of 0{Nf), and each vertex with 
the external Q and A fields is of 0{l/^/lVf). 

X)nm=i^^Q"A™ + 0(l/iVj^^), according to their power in the fields Q and A, and calculate 
all terms up to n = m = 4. An important remark is that we choose to work in a regime 
where q'^ » 1/-^/- This is consistent with our analysis in Section V, and in Appendix D, 
and still leaves the possibihty to take t ~ 1/ JNf. 



1. 0(\/iV7) terms 

The first term in Fig. 5 contributes 



SAq = - {Qr^f.{uJ = 0)+h.C.)J2 



nn&B 



X ANf 



-^4iV//?tanh(g/2T) ^ Re g„^(u; = 0), 



n 

i 



P 



^ ^ ' + 0(l/iV2) 



+ 



X 



x-4A^//3=-2/3V4A^/E^- 



(B4) 



(B5) 



where is the number of Euchdean time shces. Using the MF equation tanh(g/2T) = g/2 J, 
it is easy to show that Eqs. (B4-B5) cancel the O^^JW^) contributions in Eq. (Bl). This is 
guaranteed since the mean- field solution is a stationary point of A^'g.- 
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2. 0(1) terms 



The second diagram in Fig. 5 and the quadratic term in Eq. (Bl), give the mass terms 
of A and Q. For Q we find 



5Aq2 = 



nfieB 

1 



(B6) 



where the hnks ni/t, and 7x2/* in B' are connected by two finks in B. The form of m, (see 
also the related [37]) is 



1 1 



-m 



AT a ^ 



eioo + e) 



- = ^E 



1 



(B7) 
(B8) 



At a; = 0, and T < we get 



m = 2v — - 



J cosh^(/3g/2) 



(B9) 



which vanishes like — T near = J. In the disordered phase, where q = 0, then v = 0, 
and m2(u; = 0) = (T - Tc)/{ATTc). For ^ 0, we find that close to Tc (where g = 0) the 
difference = rn?{uj) — m^(0) is given by 



2cj 



(5 ^e^{e + uj) 13 V (e - 



E 



UJ 



E 



= ^E 



(e-a;)2e (e + a;)2e 
>0. 



Finally, the A fields are all massive 

~ E '^'^'^^+aE 772T"l2\2' 

n/xeB e It -r y ; 



(BIO) 



(Bll) 



and 5AxQ = 0. 
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3. 0{l/^/Nj) terms 



The cubic terms come from the third diagram in Fig. 5. There are contributions of the 
type Q^,Q'^X, QX^, and A^. Provided that we are interested only in the regime ~ l/^iVj, 

one needs to keep track of terms which are of 0{q^/ \jNf)- This leaves us with the following 

1 



SA; 



J2 ^(i{Qn(iMQnfiM 



"1,2 



SA: 



SA^s 



JANf niJ,eB,v 
V 

np, 

•Qn+A,i>(^i -^2)) X V^^\ 



(B12) 



(B13) 



(B14) 



where by " " we mean paths along two links, with one link in S, and by "□" we mean 

staples that begin and end on the edges of a link in B. The cubic coupling V^^^ depends on 
the external frequencies, and is given by 



e — uji 



X- 



e - a;i)2 + q"^' 



(B15) 



Next we find that SAq2x is given by 

1 



1 

X 



E 



E \Q'nfi{i^)\^{K + A„+a) 
e\u; - e) 



f3^/^ ^ {e^ + q^){{u - ey + q^) 



1 



E 



q^{e-u) 



/33/2^(e2 + 52)((^_e)2 + ^2)_ 



(B16) 



Here note that both contributions of Eq. (B16) do not involve the zero modes, since for 
a; = 0, the sums over e vanish. Next we find 



6A 



^= Re Qniti^j = 0) (A^ + A„A„+^ + A 



2 



X 



E 



qe^ 



/33/2^ (e2 + ^2)3' 



(BIT) 



and that ^^^3 — 0. 
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4. 0{1/Nf) terms 



Here we have five type of interactions terms: 6Aq4,, 6AQ3x,6AQ2x2,6AQx^,SAxi. Out of 
these, we present only 6Aq4. The reason is that 5Aq3x, and 5Aqx3 are 0{q/Nf), and that 
the 0{1/Nf) terms, 5Axi and 5Ax2q2, contribute neghgible corrections to the effective action 
of the Q fields (see Appendix C) . What we find is 



g;-(u;i+a;3-a;2)) xiyW, (B18) 



^1,2,3 



QUi^l+UJ3-UJ2))xV^^\ (B19) 



f TlfXh' 



g;-(a;i+a;3-c^2))x 21/(4), (B20) 

where " " we again mean all self-avoiding paths of length two, while the 5A^' term includes 

self-avoiding closed paths of length four, i.e. plaquettes. The quartic coupling is given by 



X 



^2 ,2 + q2 + ^^)2 + ^2 

e + uji— UJ2 



APPENDIX C: CALCULATION OF 

In this section we integrate over the massive fields Q{u} ^ and A. (The gauge transfor- 
mations are of course still massless, but do not cause any divergences in the path integral, 
since the gauge group is compact). To integrate over A we can assume that it has only a 
mass term, as close to the second term in Eq. (Bll) is negligible. The effective action for 
the link fields A^s{Q) is defined by 

exp(-^ff(g)) = exp(-^MQ.,o) X (e-^^Qx-^VA-^Vx^), (CI) 

n 

where by (, ) we mean an average with respect to SAx^-, which we evaluate by expanding the 
exponential in 1/Nf. Here we have neglected SAxi, which is of 0{1/Nf), and will contribute 

to A^siQ) only through {6Ax45Aqx2) which is of Oiq/Nf^). 

Close to Tc, only the first term in Eq. (Bll) is important, which means that only even 
terms in this expansion will survive. This leaves us with the following contributions. Prom 
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{SAqx^) we get a term of 0{q/ ^^Nf), which is hnear in 'RqQ{ijj — 0). Together with Eq. (B4) 

and the first 0{.^JJ^) term in Eq. (Bl), this yields a, 1/Nf correction for the MP equation, 
which can be absorbed in a 0{1/Nf) correction to Tc. From {5Aq2x2) wc get a con- 
tribution to the square masses of the Q fields. Since the masses of the a; 7^ fields, and of 
the zero modes are of 0(1), and 0{(f), respectively, we neglect this contribution as well. In 
fact, since the square mass of the zero modes scales like |T — T^, this contribution can also 
be considered as an 0{1/Nf) correction to Tc- Finally from {{5Aq2xY) ~ 0{1/Nf) we get 
an 0{1/Nf) contribution to the quartic interactions of the massive modes, Q((j 7^ 0). 

We now proceed to integrate out the fields Q = Q{u! ^ 0). The starting point is the 
action 

yleff(g) = E^^Q"AO - ^((5Va)')- (C2) 

n ^ 

Writing the first expression on the right hand side in terms of the zero modes — Q{oj — 
and the massive fields Q, one finds that interactions of the type ~ Oi^qj sJWf)^ and 

~ 0{1/Nf). Repeating the steps we took to integrate over A, one finds almost 
the same conclusions; there are 0{1/Nf) corrections to the square mass of the zero modes 
and to the MF equation, but there are no corrections to the quartic interactions of the zero 
modes. 

In light of the above, the effective action of the zero modes, is given by Acs{Q) 

when setting Q = 0. This gives Eq. (5.8), with the masses irinp, = m{uj = 0), and with the 
cubic and quartic terms, at zero external frequencies cui — UJ2 — uis = 0. For our discussion 
in Section V, it is sufficient to write the effective action for the zero modes above Tc. This is 
simpler, since there q — 0, which makes all the mass differences and the cubic interactions 
vanish. The result is 

1 y(^) I 1 



+2Re (^^ixK-^ix,v4>*r^-^ii+v,ix'fr^+v,i}j } ■ (C3) 

Here the second and third quartic interactions are between adjacent links lattice, and links on 
a common plaquette. Also we have = {T-T^I{ATT^, and V^^^ = ^ = 1/(48T2). 

The action Eq. (C3) has a manifest f/(l) gauge symmetry, which means that some of the 
fiuctuations will have an identically zero mass. This action is also invariant under lattice 
translations and rotations, which are the symmetries that the spin-Peierls state breaks. 
Using the notations of Fig. 3 we introduce the Fourier transform 



V^^E'^iv.e 

V -'^s JV 



1 i = 1 

gifey/2 ^ — 2 

gi(fc.+fc,)/2) ^ = 3 • (^^4) 



Here N denotes a site on the even sublattice, and k denotes a momentum in this lattice's 
first Brillouin zone. In momentum space the action looks like 
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+ X!^ii(^3 - k4)<i)k^i<i)k^i<i)k3j4>kij + Vp({fc})Re (0fcil0fe230fe320fc44) \ 
i>j J 

(C5) 

where fc belongs to the Brillouin zone of the even sublattice, and V12 = V34 = 2 cos(y(fc3 — 
fe4)/2), 1^4 = ^23 = 2cos(i;(fc3 - fc4)/2), V^ia = 2cos((i- + y){k^ - ki)/2), V^^ = 2cos((x - 
y)(fe3 - fc4)/2), and Vp = 4cos((fe2(x - y) + feay + fc4x)/2). 



APPENDIX D: CALCULATION OF THE SCALING FUNCTION f{x) 

To see the crossover behavior of Eq. (5.10) we study the susceptibihty in similar lines to 
the discussion in [10, 11]. 

/Z^0Z^0-exp(-A°,)0..0^,, 

If A^/ = 00 then ^4^^ contains only quadratic terms, and x ~ 't~^-i where t is the reduced 
temperature (T — T^/Tf.. To see deviations from this we expand the exponential in 1/Nf. 
The result is given pictorially in Fig. 6, and we find that Xij(fc, k') — X^ij^k,k', while x is 



+ 



FIG. 6: The calculation of x- The thin lines are equal to (^m^) ^, and the vertex can be read 
from Eq. (C5). 



X 



1 9 

2^ 



|m2 



X 



1 

48r24Af 



6 + 2-3)] 



1 9 

2^ 



X 



(D2) 



Here the term in the brackets comes from the vertices Eq. (C5), and also counts the number 
of fields that run in the loop. In terms of t we find that the susceptibility is 



X 



1 - 



13 1 



(D3) 



Defining x = t y ^A^/, and comparing Eq. (D3) with Eq. (5.9), we identify f{x) — 1 — 1/x 
and p — 1/2. More important, Eq. (D3) means that fiuctuations become significant when 
X ~ 0(1), indicating a crossover from MP (at x = 00) to a nontrivial fixed point of the 

effective theory in Eqs. (5.8), (5. 10). This occurs at (T - TJ/T^ ~ l/jN'f. 
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